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Recent spectroscopic observations of metal poor stars have indicated that both 7 Li and 6 Li have 
abundance plateaus with respect to the metallicity. Abundances of 7 Li are about a factor three 
lower than the primordial abundance predicted by standard big-bang nucleosynthesis (SBBN), and 
6 Li abundances are ~ 1/20 of 7 Li, whereas SBBN predicts negligible amounts of 6 Li compared to 
the detected level. These discrepancies suggest that 6 Li has another cosmological or Galactic origin 
than the SBBN. Furthermore, it would appear that 7 Li (and also 6 Li) has been depleted from its 
primordial abundance by some post-BBN processes. In this paper we study the possibility that 
the radiative decay of long-lived particles has affected the cosmological lithium abundances. We 
calculate the non-thermal nucleosynthesis associated with the radiative decay taking account both 
of the primary nuclear production reactions and the effects of secondary production as well as the 
destruction processes of energetic nuclides D, T, 3 He, 4 He, 6 Li and 7 Li. We explore the allowed 
region of the parameters specifying the properties of long-lived particles. We also impose constraints 
from observations of the CMB energy spectrum. We find that non-thermal nucleosynthesis produces 
6 Li at the level detected in metal poor halo stars (MPHSs) , when the lifetime of the unstable particles 



is of the order 



10° 



10 s and their initial abundance with respect to that of the photons is 



~ (10~ — 10~ GeV)/i? 7 o, where _B 7 o is the emitted photon energy in the radiative decay. We 
conclude that the most likely nucleosynthetic scenario involves two different processes. First, a 
non-thermal cosmological nucleosynthesis associated with the radiative decay of unstable particles; 
and second, the stellar depletion of both of the primordial lithium isotopic abundances. 



I. INTRODUCTION 

In standard cosmology, the universe is thought to have 
experienced big-bang nucleosynthesis (BBN) at a very 
early stage. D, T, 3 He, 4 He, 6 Li, 7 Li and 7 Be are pro- 
duced appreciably at this epoch. Different types of ob- 
servations have been made for the purpose of determin- 
ing the primordial elemental abundances. These observa- 
tional results provide rich information about the cosmic 
chemical evolution of the light nuclear species. Today 
BBN provides a very precise tool for inferring the con- 
dition of the early universe. BBN in the standard cos- 
mology explains relatively well the inferred primordial 
abundances for a narrow range of the universal baryon- 
to-photon ratio rj. 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite has measured the temperature fluctuations of the 
cosmic microwave background (CMB) radiation, and pa- 
rameters characterizing the standard big bang cosmology 
have been deduced [l|,|2| f rom these data. For the baryon- 
to-photon ratio r\ deduced from fits to the CMB, the BBN 
model predicts abundances of the light elements which 
are more-or-less consistent with those inferred from as- 
tronomical observations. This near agreement places sig- 
nificant limits on non-standard models which influence 
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the cosmic nuclear abundances. 

In this regard, unstable massive particles decaying dur- 
ing or after the BBN epoch are strongly constrained 0, 
QLIEE1' For example, radiative decay pi El 01 ' hadronic 
decay, or annihilation 0, 0, 0] have been studied re- 
cently, and several critical constraints on the relic par- 
ticle properties were derived. These particle processes 
induce electromagnetic and/or hadronic showers which 
lead to the destruction of preexisting nuclei and to the 
production of different nuclear species. In turn, these 
modifications to the light element abundances are used 
to constrain theories for the decay of relic particles. The 
physics associated with the general decay process (includ- 
ing a hadronic decay component) is explained in [T^.ll3||. 

Super-weakly interacting massive particles (SWIMPs) 
have been proposed 01 as candidates for cosmological 
dark matter. They are a form of non-baryonic dark 
matter derived from theories beyond the standard model 
(SM) such as supersymmetry or universal extra dimen- 
sions (UED) [15| . The stable particles in these theories 
may constitute the dark matter, while the unstable par- 
ticles could decay to modify the nuclear abundances in 
the early universe. This type of dark-matter particle in- 
teracts super-weakly with the SM fields and cannot be 
detected directly in experiments searching for conven- 
tional WIMPs. Thus, further studies of decay-induced 
nucleosynthesis could help to constrain such new particle 
theories. 

Spectroscopic lithium abundances have been detected 
in the atmospheres of metal poor stars. Nearly constant 
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abundances of 6 Li and 7 Li in metal-poor Pop ulation II 
(Pop II) stars have been inferred [lj| [lTJ . Spectro- 
scopic measurements obtained with high resolution in- 
dicate that metal poor halo stars (MPHSs) have a very 
large abundance of 6 Li, i.e. at a level of about a twenti- 
eth that of 7 Li. This is about three orders of magnitude 
larger than the SBBN prediction of the 6 Li abundance. 

It was suggested some time ago p| that a high 6 Li 
abundance could be produced via particle decay. More 
recently, radiative particle decay in particular has been 
proposed Q as the source of the abundant 6 Li observed in 
low metallicity stars. Furthermore, it was suggested 
that hadronic decays which occur at around 1000 seconds 
after the big-bang could simultaneously resolve both the 
6 Li and 7 Li abundance issues. This probability has been 
studied within the context of supersymmetry, and there 
is sufficient parameter space for this to be possible . A 
previous study on the radiative decay, however, was un- 
able to provide a thorough solution to the lithium abun- 
dance puzzle |2(j. In this paper, therefore, we consider 
this issue to be as of yet unresolved. We independently 
calculate the nucleosynthesis triggered by the radiative 
decay processes of long-lived relic particles over a wide 
range of parameters specifying the properties of the relic 
particles. We take into account the primary, secondary, 
and tertiary processes resulting from the electromagnetic 
cascade showers which both produce and destroy the light 
elements. We then constrain the abundance of long-lived 
particles from the calculated nucleosynthesis. We do not 
find, however, a simultaneous solution to both the 7 Li 
and 6 Li abundances unless there is stellar destruction of 
lithium. 

This paper consists of the following structure. In 
Sec. [n] we describe briefly the formulae of non-thermal 
nucleosynthesis including the relevant photo-dissociation 
and nucleus-nucleus collisions. In Sec. IIIII we summa- 
rize the observed light element abundances and the con- 
straints adopted in this work. The calculated results are 
shown in Sec. II VI where we derive the constraints on the 
parameters of decaying particles. In Sec. the con- 
straints from the CMB spectrum are imposed, and the 
effects on small scale structure are analyzed. The most 
interesting parameter region is deduced by taking the ob- 
servational constraints on the light element abundances 
into account. We conclude that our model can explain 
the desired 6 Li production by non-thermal nucleosynthe- 
sis even if there is stellar destruction of both lithium iso- 
topes to explain the observed 7 Li. 



II. MODEL 

We assume the creation of high energy photons from 
the radiative decay of a massive particle with a lifetime of 
10 2 - 10 12 s. In this section we describe equations to cal- 
culate the non-thermal nucleosynthesis triggered by the 
appearances of high energy photons from decay. See 
for details on the formulation which we adopt. 



A. Initial non-thermal photon spectrum 

We assume that the decaying dark particle is non- 
relativistic, and almost at rest in the expanding universe. 
We denote the imaginary particle by X, [e.g. corre- 
sponding to the next-to-lightest supersymmetric particle 
(NLSP)] with a mass Mx that decays into a photon plus 
another dark-matter particle, [e.g. the lightest supersym- 
metric particle (LSP)]. We represent the emitted photon 
energy by E l0 . 

When an energetic photon emerges, it interacts with 
the cosmic background and induces an electromagnetic 
cascade shower. The faster processes are pair production 
through background photons 7b g (77bg — * e + e~) and 
inverse Compton scattering of produced electrons and 



positrons through background photons (e 7b g 



=7). 



These two processes produce electromagnetic showers 
and the non-thermal photon spectrum realizes a quasi- 
static equilibrium 0, [2^ • The attained zeroth genera- 
tion photon spectrum can be written p^ |. 



Pi{ E i) 



KoiEx/E^) 1 - 5 for E 1 < E x 
K (Ex/E 7 ) 2 for Ex < E^ < E c 



(1) 







for E c < E. 



where K = E 7 q/ (E x [2 + ln(Ec/Ex )]) is a normaliza- 
tion constant fixed by energy conservation of the injected 
photon energy. This spectrum has a break in the power 
law at E 1 = Ex and an upper cutoff at E 1 = Ec- We 
take the same energy scaling with the temperature T of 
the background photons as in |2l|, i.e. Ex = m 2 /80T 
and Ec = m 2 /22T, where the cascade spectrum was cal- 
culated by numerically solving a set of Boltzmann equa- 
tions. 

Above Ec = m 2 /22T, the rapid interaction between 
photons and electrons or positrons causes almost all en- 
ergetic photons to quickly lose their energy. This leaves 
effectively no number spectrum for energies E 7 > Ec- 
In the context of the theory of electromagnetic cascade 
showers j2^|, the scale Ex — m 2 /80T corresponds to the 
point where the energy of the energetic photon is below 
the threshold for double photon pair creation. Hence, 
pair creation from interactions between energetic pho- 
tons and the cosmic background radiation (CBR) ceases 
to operate. 

The zeroth-generation non-thermal photons experi- 
ence additional processes including: Compton scattering 
(7 e be ~~ * 7 e± )i Bethe-Heitler ordinary pair creation in 



nuclei (jN^g 



~e N); and double photon scattering 



(77bg —* 77)- These slower processes further degrade the 
quasi-static equilibrium photon spectrum. 

Because the rates of these electromagnetic interactions 
are faster than the cosmic expansion rate, the photon 
spectrum is modified into a new quasi-static equilibrium 
(QSE). This distribution is given by 



A/? SE (£ 7 



n x p 1 {E 1 ) 
T y (Ej)Tx ' 



(2) 
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where nx = n x (l + z) 3 exp(— t/rx) is the number den- 
sity of the decaying particles at a redshift z, and tx is 
its mean life. The quantity T 7 is the energy degradation 
rate through the three slower processes for the zeroth- 
generation photons. We use this steady state approxima- 
tion for the cosmic non-thermal constituent of photons. 
The Boltzmann equation describing the electromagnetic 
cascade processes is compiled in [21|. We used their rel- 
evant reaction rates to compute the cascade shower. 

Because the cutoff scale of the photon energy is in- 
versely proportional to the background temperature T, 
the cutoff scale comes to large as the universe cools. As a 
result, the number of the energetic non-thermal photons 
increases at low T. 



where we have defined the reaction rate 
n°Cx / 1 \ 3/2 

x f dE, (-^Af^iE,)] a.+B^ciE,), 
Jo \E l0 nx J 

(5) 

where C,x = rE-yo. 

The primary reactions and their cross sections are 
taken from |jj . 



C. Secondary nucleosynthesis 



B. Primary nucleosynthesis 

The equation for the production and destruction of 
nuclei by non-thermal photons is given by 



dY A _^ N A (T) 

T 



poo 



N A (P) 
N A (P)\ 



/ dE^ SE (£ 7 ) ct 7+ (£ 7 ) , 
Jo 



(3) 



where Yi = ni/np is the mole fraction of a particular 
nuclear species i, and rii and np are number densities 
of nuclei i and baryons. The first and second term on 
the right-hand side are the source and sink terms for 
nucleus A. The source term includes the creation of nu- 
clide A from all possible target nuclides T through a non- 
thermal photo-dissociation reaction 7 + T — > A. And the 
sink term includes the destruction of nuclide A that is 
specified by the reaction 7 + A — > P for any produced 
nuclides P. E~ is a non-thermal photon energy. The 
cross sections of processes 7 + T — > A and 7 + A — > P 
are denoted by a 7 +T~tA{E 7 ) and a,+A^p(E- ( ). Further 
we use Nk(L) to represent the number of a particular 
particle K relevant to a process 7 + a — > b, where K is 
either a or 6 and L is T or P. For example, in the pro- 
cess 4 He(7,d)c?, iVd( 4 He)=2 because two deuterons are 
produced from 4 He. 

Next, we formulate the rate equation. We define 
r = n x /n® and H r = ^SirGp^/S, where the super- 
script denotes present values (z = 0), therefore n° 
and /9° ad are the present CBR photon number density 
and present radiation energy density, respectively. Then 
nx = n x (l + z) 3 exp(— t/rx) is transformed into nx = 
n x (l/(2H r t)) 3 / 2 exp(— t/rx), so that the rate equation 
© becomes 



Y A 



N A (py 



[Mi 



Yp 



N P (P)\ 



(4) 



If the photo-dissociated light nucleus of a primary re- 
action has enough energy to induce further nuclear reac- 
tions, then secondary or tertiary processes are possible. 

The equation describing the secondary production and 
destruction is obtained by taking into account the en- 
ergy loss of nuclear species while propagating through 
the background. In most situations, the energy loss rate 
is faster than any sink of primary nuclei so that they 
can be ignored in the evolution of the primary particles. 
However, for unstable particles, the sink terms must be 
evaluated. In general, because of the high reaction rate, 
the primary particles establish a quasi-static equilibrium. 
The abundance evolution is then represented by 



dY s 
dt 



N A (A) N s {S) 
N t {A)\ N a {S)\ 



Y T Y T >[TT'] S - (sink term), 



T.T' 



( 6 ) 

where N^{L) has the same meaning as above. The quan- 
tity [TT']s is given by, 



[TT'} 



v(<fu ( 1 
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<7a+t'->s(Ea)Pa 







b A (E A ) 
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TX 



x exp 
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£a(E-,) 
Ea 



dE, 



Af^ SE {E y ))a 7+ T^A(E 1 
^a(Ea) 



b A {E; 



(7) 



This is the reaction rate for a secondary process 
T(7, Xi)A(T', X2)S with any combination of particles 
Xi, A, and X2- The quantity r\ is the baryon-to photon 
ratio: n = n B /n®, and Pa is the velocity of the primary 
particle A. bA = —dE/dt is the rate of energy loss of 
the primary particle. The energy loss rates are taken 
from 01 • Coulomb scattering (Ne^ s — > Ne ± ) is the 
dominant loss process for primary nuclei AT. Then Ta 
is the coefficient of the primary particle sink, and in the 
case of an unstable nuclide, its value is the decay rate. 
£a(Ej) is the energy of the nuclide A produced by the 
photo-dissociation process 7 + T — > A, simultaneously 
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£ a ~ 1 (Ea) is the energy of the non-thermal photons which 
produce the primary species A with energy E A . 
To simplify Eq. Q we define 



TX 



E^qUx 



A/? SE (£ 7 ), 



and 



P(E A ,E^=exp 



£a{E-,) 



dE 



E A 



ta(e; 

b A (E A 



(8) 



(9) 



Then we denote the factors preceding the first integral 
by an overall normalization B to obtain 

[TT%=B f A{E iE A ° A ^ EA ^ A 

Ec 



b A {E A ) 

JS-^Ea) 

(10) 



where lower limit for the external integral is replaced with 
the threshold energy of primary nuclide for the secondary 
reaction E p ^h- Since the photon spectrum includes few 
energetic photons above Ec, the upper limit for internal 
integral can be changed to Ec- £a{E 1 ) and £^ 1 {Ea) 
are derivable in the limit of low energy scattering, where 
the relevant nuclei are non-relativistic. The adoption of 
this limit is reasonable because most of the reactions oc- 
cur with low energy photons near threshold. Thus, the 
produced primary nuclei are mostly non-relativistic. 

If the primary product particles are further destroyed 
by interactions with background particles, this secondary 
destruction process can have a large contribution to non- 
thermal nucleosynthesis. Thus we also take into account 
the destruction of D, T, 3 He and 6 Li after primary pro- 
duction by abundant background nuclides. For the de- 
struction processes d+p and 3 He+p, we used the cross 
sections from 0, |2^. For the t(p,dp)n and t(p,p2n)p 
reactions, we used cross sections from the mirror nucleus 
reactions 3 He(p,dp)p and 3 He(p,2pn)p, respectively. Al- 
though the cross sections of t+p and 3 He+p reactions 
are different because of the Coulomb interaction, these 
mirror nucleus reaction cross sections are almost identi- 
cal l2al. We adopt the cross section for 6 Li(p, 3 He) 4 He 
from |27| . 

It has been found that only 6 Li can be produced by 
secondary non-thermal nucleosynthesis at a much higher 
level than by the SBBN, while the production of other 
nuclei e.g. D, 7 Li and 7 Be are insignificant 0, 0, IT^ |. 
The relevant processes in the secondary non-thermal pro- 
duction of 6 Li involve interactions of background 4 He 
with primary tritium or 3 He particles. Threshold en- 
ergies of the 4 He(t,n) 6 Li and 4 He( 3 He,p) 6 Li reactions 
are 8.3870 MeV and 7.0477 MeV, respectively. We have 
taken into account these two reactions with their cross 
sections from |9j|. 



The survival probability should also be considered in 
order to calculate the precise abundance I12J. This de- 
struction of 6 Li is included in the SBBN code There- 
fore the formulation of the survival probability should 
not include the process after the thermalization of 6 Li to 
prevent a double counting of 6 Li destruction. 



III. OBSERVATIONS OF LIGHT ELEMENT 
ABUNDANCES 

A. Light element abundances 

The primordial abundances of D, 3 He, 4 He, and 7 Li 
are inferred from various observations. Recently 6 Li has 
been measured for numerous old halo stars and a so- 
called 6 Li plateau as a function of metallicity appears to 
exist pflj . The reliability of inferred primordial abun- 
dances from several types of observations is difficult to 
estimate because of systematic errors. Here, we summa- 
rize the observational data and our adopted constraints. 

Deuterium is measured in absorption spectra at high 
redshift toward QSOs. Since the ratio D/H is low, D 
can only be detected in absorption systems with high HI 
column densities. This causes a difficulty for such obser- 
vations. From the Keck I HIRES spectra of Q1243+3047, 
a deuterium abundance of D/H= 2.42l°;|g x 10 -5 is es- 
timated [23 . The best estimate of the primordial D/H 
from absorption systems toward five QSOs reported in 
that paper, is 



D/H = 2.78: 



u.:SH 



x 10" 



(11) 



However, the five values have a larger dispersion than 
that inferred from the individual measurements. Hence, 
the errors are probably underestimated. We therefore 
take the highest value of D/H in the five values, D/H= 



3.98^ 6? x 10~ 5 and allow for a two sigma increase above 
this value to fix an upper limit. In the standard theory 
of galactic chemical evolution, deuterium continually de- 
creases in time by the processing of interstellar material 
through stars. Hence, we take the present abundance 
of deuterium as the lower limit. Column densities of DI 
along seven lines of sight have been estimated from obser- 
vations with the Far Ultraviolet Spectroscopic Explorer 
(FUSE), and the local interstellar medium at a distance 
of 37-179 pc has been probed [3(j ■ The weighted mean 
value of DI/HI for five data, with reliable values for col- 
umn densities of HI determined, is (1.52 ± 0.08) x 10~ 5 . 
When we allow for a two sigma lower limit, we get the 
following constraint on the primordial abundance of deu- 
terium 



1.4 x 10~ 5 < D/H < 5.2 x 10" 



(12) 



3 He is measured in Galactic HII regions by the 
8.665 GHz (3.46 cm) hyperfine transition of 3 He + [3l| . 
A plateau with a relatively large dispersion with respect 
to metallicity has been found at a level of 3 He/H=(1.9± 
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0.6) xlO -5 . There is a problem, however. Although stars 
are thought to produce 3 He, and 3 He enhancement is ob- 
served in planetary nebulae, the chemical evolution of 
3 He has not been detected in the Galaxy during the last 
4.5 Gyr [32,133]. This fact has recently been confirmed by 
more precise determination of the helium isotopic com- 
position of the local interstellar cloud (34[. It is not yet 
understood, therefore, whether 3 He has increased or de- 
creased through the course of stellar and galactic chem- 
ical evolution (3^, |3|J. Thus, we adopt the two sigma 
upper limit from Galactic HII region abundances, that is 

3 He/H < 3.1 x 10~ 5 . (13) 

We will discuss later an implication for a tighter con- 
straint as suggested in [34|, i.e. 3 He/H< (1.6 ± 0.3) x 
10~ 5 , too. We do not give a lower limit due to the large 
uncertainty in the galactic production of 3 He. 

4 He is measured in the HII regions of metal-poor ex- 
ternal galaxies where chemical evolution is thought to 
be minimal. The primordial abundance is estimated to 
be Y = 0.2421 ± 0.0021 by extrapolating the abundance 
to zero metallicity (O/H=0) 0. However, in \^ it is 
noted that there are sources of systematic uncertainty 
in determinations of the 4 He abundance. They sug- 
gest somewhat larger error bars with an abundance of 
Y = 0.249 ± 0.009. They thus adopt a primordial 4 He 
abundance within the conservative range of 

0.232 < Y < 0.258. (14) 

We also adopt this constraint for the primordial 4 He 
abundance. 

7 Li is measured in metal-poor halo stars (MPHSs) by 
the spectra of their atmosphere. There is about a factor 
of three under-abundance of 7 Li in MPHSs with respect 
to the SBBN prediction when using the baryon-to-photon 
ratio 77 inferred from the analysis of the CMB anisotropy. 
This is called the lithium problem 0, llTllpl . For a re- 
cent review of the lithium problem, see |39| . where pos- 
sible resolutions of this problem are discussed. 

Recently, high quality spectra for 24 metal-poor halo 
dwarfs and sub-giants have been obtained j2^] . They es- 
timated a mean 7 Li abundance of \0ge7u — 2.21 ± 0.07. 
Different groups derive somewhat different values. We 
adopt the estimate of [23 for the primordial 7 Li abun- 
dance allowing for a two sigma range, loge7Li = 2.07 — 
2.35 or 

1.17 x 10~ 10 < 7 Li/H < 2.23 x 10~ 10 . (15) 

In addition, we take into account the possibility of the 
stellar depletion of lithium up to 0.5 dex to derive an 
upper limit to the primordial abundance, 

1.1 x 10~ 10 < 7 Li/H < 7.1 x 10~ 10 . (16) 

As for the depletion of lithium in halo stars, it has 
been reported |40j that the mean lithium abundance and 
its dispersion appear to be lower for dwarf stars than 



for turn-off and sub-giants. 7 Li abundances of 28 such 
halo sub-giants have been measured j^l). The result 
is that (excluding the extremely lithium-rich sub-giant 
BD +23 3912) the maximum abundance is loger Li = 
2.35. This is well below the SBBN predicted value. A 
calculation of the stellar depletion of lithium isotopes by 
atomic and turbulent diffusion leads [421 ] to an abun- 
dance reduction by a factor of at least 1.6-2.0 of the 
7 Li abundance for Population II stars with metallicity 
[Fe/H] < -1.5. 

6 Li has also been measured in MPHSs by spectroscopy. 
In [23, 6 Li was detected at a better than two sigma sig- 
nificance in nine of the 24 stars observed. They suggest 
that a 6 Li plateau exists at loge6 Li w 0.8. This plateau 
implies important information because the SBBN would 
predict that much less primordial 6 Li be formed ( 6 Li/ 7 Li 
~ 10~ 5 ). Therefore, some mechanism should have pro- 
duced almost all 6 Li in MPHSs. There are several candi- 
dates. For example, gravitational shocks during the hi- 
erarchical structure formation of the Galaxy could have 
accelerated cosmic rays (CRs) and produced 6 Li by a — a 
fusion reactions at a very early stage |43l 13 ] . A scenario 
has also been suggested 0, 0] whereby cosmological 
CRs (perhaps related to pop III stars) were produced 
at an extremely early epoch up to the formation of the 
Galaxy. These CRs could then have produced 6 Li by 
a + a fusion. See the discussion in |46] for a recent sum- 
mary of various candidates for 6 Li production. Another 
new mechanism has been proposed recently |47j whereby 
6 Li originated from the binding of negatively charged par- 
ticles to background nuclei. Since multiple processes have 
possibly synthesized 6 Li at an early epoch, we do not put 
limits on the primordial abundance of 6 Li. However, we 
adopt the average value of the abundance derived from 
the eight MPHSs with detections as a guide to constrain 
possible production from radiative decay, 

6 Li/H w 6.6 x 10~ 12 . (17) 

B. Cosmic microwave background anisotropy 

Very precise data have been obtained by observa- 
tions of the spectrum of temperature fluctuations in the 
CMB. The WMAP data have been analyzed and the en- 
ergy density of baryons in the universe have been de- 
duced along with other cosmological parameters. This 
leads to VL b h 2 = 0.0224 ± 0.0009 for the WMAP first 
year data Q and tt b h 2 = 0.0207^;^ for the WMAP 
three year data in the running scalar spectral index 
model. We adopt the first year result, corresponding to 
V = (6-l±o;|) x 10~ 10 . The SBBN with the WMAP fl b h 2 
parameter region has been calculated including the un- 
certainties of the inferred fl hh 2 and the the reaction rate 
uncertainties on the SBBN |4||. Their result is: 

D/H = (2.60tgi?) x 10_5 ( 18 ) 
3 He/H = (1.04 ± 0.04) x 10" 5 (19) 
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Y = 0.2479 ± 0.0004 (20) 
7 Li/H = (4.15±g;t|) x 10- 10 . (21) 

IV. RESULTS 

We have calculated the cosmological nucleosynthesis 
including the SBBN and non-thermal nucleosynthesis in- 
duced by the radiative decay of a long-lived massive 
particle. The SBBN was computed using the Kawano 
code |49j with the use of the new world average of the 
neutron lifetime We added the non-thermal com- 

ponents of Eq. 10} and © also taking into account the 
destruction of secondarily produced 6 Li. We deduced 
properties of the decaying particle in terms of its lifetime 
Tx and abundance with the baryon-to-photon ratio 
fixed at X] = 6.1 X 10~ 10 . We checked the effect of sec- 
ondary destruction of the primary non-thermal nuclides. 
We confirmed that the energy loss is much faster than 
the destruction (see the discussion for non-thermal 6 Li 
production of 0)- The secondary destruction processes 
of primary nuclides were not very efficient (destruction 
probabilities are < 0(1O -3 )). The final abundances were 
therefore not much affected by secondary destruction. 

The point is that the time scale of the Coulomb loss for 
the non-thermal nuclides is much smaller than those of 
the destruction reactions at low energies. Furthermore, 
the destruction probabilities of non-thermally produced 
6 Li are very small, and produced 6 Li survives |12j . 

We have derived the constraints on the lifetime tx and 
abundance parameter C,x from the adopted limits for the 
cosmological light element abundances. Our result is very 
similar to that of 0, since we use the same formulation 
for non-thermal nucleosynthesis and adopt their fitted 
cross sections. A detailed explanation has been given 
in 9] for the systematics of the radiative decay. 

Fig. n shows a contour of the 4 He mass fraction Y > 
0.232 (red line) in the (tx,(x) plane. Above this con- 
tour, Y < 0.232. Since 7 Li is more weakly bound than 
4 Hc, the destruction happens even in shorter lifetime 
conditions. The contour of 7 Li for the lower limit of 
7 Li/H> 1.1 x 10~ 10 (blue line) has the shape plotted 
on Fig. ^ Contours for D/H upper and lower limits, 
D/H< 5.2 x 10~ 5 (green solid lines) and D/H> 1.4 x 10~ 5 
(green dashed lines), respectively are also shown. Con- 
tours for the 3 He/H upper limit, 3 He/H< 3.1 x 10~ 5 
(black lines) are also drawn. The shape of this contour 
is very analogous to that for the D/H upper limit. This 
reflects the fact that the photodisintegration of 4 He is the 
main cause of the production of both D and 3 Hc. Finally, 
the contour for the MPHSs level of 6 Li/H=6.6 x 10~ 12 is 
plotted to see the behavior of the 6 Li non-thermal pro- 
duction (orange line). 

Fig. El shows the derived constraint on tx and (x for 
an unstable particle from the consideration of the light 
element abundances above described in a model with 
r\ = 6.1 x 10~ 10 . All limits from the light element abun- 
dances used in this study, for D/H, 3 He/H, Y and 7 Li/H 




He over 



t x (sec) 

FIG. 1: Contours in the (tx , (,x) plane corresponding to the 
adopted constraints for the primordial abundances in mod- 
els with a fixed value of r\ = 6.1 x 1CU 10 . Contours for the 
mass fraction of 4 He Y — 0.232 (red line) and the number 
ratios of 3 He/H=3.1 x 1(T 5 (black lines), D/H=5.2 x 1CT 5 
(green solid lines), D/H=1.4 x 10 -5 (green dashed lines), and 
7 Li/H=l.l x 10 _ 10 (blue line) are shown. The contour of 
6 Li/H=6.6 x 10 -12 (orange line) is also drawn. The notation 
"over" and "low" identifies overproduced and underproduced 
regions, respectively. 




t x (sec) 

FIG. 2: Colored regions identify the excluded area in the 
parameter space (tx,Cx) for models with a fixed value of 
ri = 6.1 X 10" 10 . The dark shaded region is excluded by 
an overabundance of 3 He, whereas the light shaded region is 
mostly excluded by an underabundance of deuterium. 

have been taken into account. The 3 He overabundant 
region is shaded by the dark color, and the rest of the 
excluded region a light color. The light colored region 
is fixed largely by the deuterium underproduction. For 
Tx^lO 6 s, 3 He provides the strongest limit on the abun- 
dance parameter yielding, 

Cx < 4 x 10~ 13 GeV , (22) 

while for shorter lifetimes (tx ~ 10 4 — 10 6 s) the limits 
are from D implying (x ^10 -9 GeV. 

Uncertainties from the reaction cross sections have 




been discussed in |9j. The uncertainties from the pho- 
todissociation and secondary cross sections are of the or- 
der of 10 %. In addition, we confirmed that effects of 
errors in the cross sections for the secondary destruction 
reactions are negligible. The BBN light clement abun- 
dances are then the dominant source of uncertainty. The 
BBN cross section uncertainties are not included in our 
study, but they are estimated [e.g. Eq. I|18|)- 1)21(1] . and 
one can apply those to our non-thermal calculation re- 
sults. 



V. DISCUSSION 

A. Distortion of the CMB spectrum 

The abundance of massive decaying particles is also 
constrained from the observed spectrum of CMB radia- 
tion. Since non-thermal photons produced by the radia- 
tive decay deform the blackbody spectrum of the CMB, 
this is limited by the consistency of the observed CMB 
data with a blackbody spectrum 0, For epochs 

earlier than z ~ 10 7 , thermal bremsstrahlung, [i.e. free- 
free emission (eN — > eNj), where N is an ion] and 
radiative-Compton scattering (e~7 — > e~"/"/) act effec- 
tively to erase any distortion of the CBR spectrum from 
a blackbody. For the decay in epochs 10 5 < z < 10 7 , 
processes changing the photon number become ineffec- 
tive, and Compton scattering {^e~ — > 7e~) causes the 
photons and electrons to achieve statistical equilibrium, 
but not thermodynamic equilibrium. Then, the photon 
spectrum should have a Bose-Einstein distribution 



1 



e e 7 /T+jU _ I 



(23) 



where ji is the usual dimcnsionless chemical potential de- 
rived from the conservation of photon number. 

Analyses of the CMB data suggest a relatively low 
baryon density so that double Compton scattering domi- 
nates the thermalization process. For small energy injec- 
tion from the radiative decay, the chemical potential can 
be approximated analytically [fill IH^ as 



(j, = 4.0 x 10" 



where 



T dc =6.1xlO b s 



4 




1/2 


f Cx 1 




-10 6 s- 




10" 9 GeV 



To 


-12/5 


'n b h 2 ' 


4/5 


'i -y/2" 


4/5 


2.725 K 




0.022 




0.88 





(24) 



(25) 



where To is the present CMB temperature, and h = 
ifo/(10G km s _1 Mpc -1 ) is the normalized Hubble pa- 
rameter. 

For a late energy injection at z < 10 5 , Compton scat- 
tering produces little effect and cannot establish a Bose- 
Einstein spectrum. The distorted spectrum is then de- 
scribed by the Compton y parameter. There is a rela- 
tion between y and the amount of the injected energy, 




100 1000 10* 10 J 10" 10' 10° 10° 10 1 " 10 11 10" 
t x (sec) 

FIG. 3: Same as in Fig.|5]except that the black shaded region 
is superimposed. This region shows the region excluded by 
the consistency requirement of the CMB with a blackbody. 
The white band identifies parameter values where the free- 
streaming of decay products would lead to a suppression of 
small scale structure growth in a model with Qyh 2 ~ 0.11. 
The curved line identifies the contour of 6 Li/H = 6.6 x 10 -12 , 
corresponding to the abundance of 6 Li observed in MPHSs. 
The region above the contour and below the nucleosynthesis 
and CMB constraints is allowed and abundant in 6 Li. The 
contour for 6 Li enhanced by three times as much MPHSs value 
is also drawn (dashed line). The solid line corresponds to the 
sensitivity of the planned ARCADE mission as labeled. 



AE/Eqbr = 4y, where AE and -Ecbr are the total en- 
ergy injected and the CBR energy, respectively. The ra- 
tio of the energy injected by radiative decay to the CBR 
energy per comoving volume can be expressed as 



E 



AE 



Ecbr 2.7T(i c ff) 



n-xo 



n 7 o 



c 



A" 



2.7T(t cff ) 



(26) 



where t cS = [r(l + /3)] 1//3 t x for T oc i"' 3 , with T(x) the 
gamma function of argument x. The factor of 2.7 comes 
from the fact that the average energy of the CBR photons 
at a given temperature is 2.7T. The time-temperature 
relation is T oc t~ x / 2 for a radiation dominated universe. 

The CMB spectrum has been well measured and the 
deduced limits are < 9 x 10" 5 , \y\ < 1.2 x 10" 5 jH 
and n b h 2 - 0.022 with h - 0.71 pj. Therefore, the high 
abundance parameter region of £x is excluded by the /i 
and y limits. In Fig.[3]the black shading indicates the pa- 
rameter region excluded by the CBR distortion limit. For 
a lifetime shorter than tx — 4 x 10 11 s flbh 2 ~ 8.8 x 10 9 s, 
the decay is constrained by the chemical potential /i. On 
the other hand, when an unstable particle decays later, 
the CBR spectrum is limited by the Compton y param- 
eter. The parameter region of relatively long lifetime 
(10 10 s < tx) is found to be constrained by the CMB 
spectrum more strongly than the light element abun- 
dances. 

The discussion so far has assumed that rj at the BBN 
epoch can be determined from the power spectrum of 
CMB temperature fluctuations. However, r\ could have 
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changed from the BBN epoch to the time of the recom- 
bination |54| . The entropy production resulting from the 
photon emission by radiative decay has been estimated 
and it is very small for a large part of the region al- 
lowed by the non-thermal nucleosynthesis and CMB dis- 
tortion constraints in Fig. |21 However, for very early 
decay (tx < 10 4 s) entropy production could be large. 
An entropy change by as large as AS/Si ~ 1CP 1 seems 
unlikely since the CMB-favored value of 77 = 6.1 x 10~ 10 
is transformed into 6.7 x 10 -10 at the BBN epoch and this 
would bring about a more severe discrepancy of the 7 Li 
abundances between the SBBN prediction including the 
non-thermal processes and estimations from the MPHSs 
observations. See Fig. 3 in [54| 



B. Effects on small-scale structure 

The decay products obtain momentum from the decay 
by simple momentum conservation. This momentum can 
affect structure formation in the universe. In this regard 
there are useful constraints on the present root-mean- 
square velocity of warm dark matter from observations 
of the Lyman-a forest |55| and super-massive black holes 
at high redshift [5||. Limits from the early reionization 
have also been reported [5(| |13 • The free-streaming de- 
cay products must not cause an inconsistency with ob- 
servations. This consideration of free-streaming has been 
discussed in detail in |igj|. 

Cold, collision-less, non-self-interacting dark matter 
explains the trend of the observed structures larger than 
~ 1 Mpc. However, it does not successfully explain the 
structure at smaller scales. As a solution to the small 
scale structure problem, the decay of SWIMPs has been 
proposed |58|. The velocity dispersion reduces the 6- 
dimensional phase space density and prevents the forma- 
tions of cuspy halos. Furthermore, the free-streaming of 
dark matter particles damps the power spectrum at small 
scales. They discuss the paradigm in which the dark mat- 
ter is composed mostly of the decay products. In this case 
the parameters which lead to reasonable free-streaming 
predict roughly the correct phase space density. There- 
fore we only discuss the free-streaming scale. 

If the free-streaming scale is reasonably small ~ 1 Mpc, 
the free-streaming scale Afs from the decay time to the 
present (z ~ 0) is described [5gj as 



Afs ~ 1-5 Mpc uy 



TX 



10 6 



1/2 



(27) 



where uy = py /m is momentum at the decay time di- 
vided by the mass of the particle Y produced by the 
decay. In the range of tx < 10 12 s, this approximation 
involves a maximum error of only a factor 1.5 in the nu- 
merical factor. The largest error comes when tx = 10 12 s 
and correspondingly uy ~ 10~ 3 . From momentum con- 
servation during the decay, the momentum of product py 
is equal to the emitted photon energy E^q. The param- 



eter £x is written as 

(x = 2.6 x 1CT 8 GeV (n Y h 2 )u Y , (28) 

where fly is the energy density normalized to the critical 
density. 

Eqs. 1)2 7|l and l|28|l connect the tx and (x by the elim- 
ination of uy. Such decay can be a resolution to the 
small scale structure problems if the free-streaming scale 
is 0.4 Mpc < A FS ^ 1-0 Mpc [H^j. This limit determines 
the suitable parameter space for tx and Cx under the 
assumption of Qyti 2 ~ 0.11 0. This region is shown in 
Fig.|31as a white band bounded by solid lines. The region 
above this band is excluded in this case, since the decay 
products erase the larger scale fluctuations. The region 
below this band will be discussed in the next subsection 
in connection with the 6 Li abundance. Apparently, the 
decay at tx ^ 10 5 s is not realistic because the effect on 
the light element abundances or CBR spectrum is pro- 
hibitively large. On the other hand relatively early decay 
at tx < 10 5 s is viable. Applying this result to a model 
for a gravitino SWIMP with a photino NLSP [5^|, we 
obtain toswimp < 200 GeV and Am ~ 400-700 GeV. 



C. 6 Li-producing parameter region 

In this paper we analyze the possibility that the ra- 
diative decay of long-lived particles produces 6 Li by non- 
thermal process while having almost no effect on 7 Li or 
other nuclides produced in the SBBN. Ellis, Olive & Van- 
gioni researched the possibility that the radiative decay of 
unstable particles explains the discrepancy of the BBN 
calculated 7 Li abundance and low 7 Li plateau derived 
from observations [2C| . They found that in the param- 
eter region where 7 Li is photo-dissociated down to the 
level of the 7 Li plateau, either the D abundance was too 
low or the ratio 3 He/D was unacceptably large in the con- 
text of standard stellar evolution and chemical evolution. 
Consequently, they concluded that radiative particle de- 
cays cannot be a cause for the 7 Li abundance difference. 
They also mentioned the possibility of 6 Li production in 
their paper. 

Uncertainties remain in estimations of the Li abun- 
dance in stellar atmospheres, and the probability of de- 
pletion in stars has not been excluded. The difference 
between abundances determined by different analysis ap- 
proaches is also a concern. Therefore, we suppose that 
the discrepancy of the 7 Li abundance is caused by stellar 
depletion or some other systematic effect. Then the 6 Li 
abundance in the early universe should have been larger 
when first engulfed in a star than the value presently de- 
duced from observations of MPHSs. Assuming that is 
the case, we impose the following constraint on the 6 Li 
abundance after the radiative decay process, 



3 Li/H > 6.6 x 10' 



-12 



(29) 



If this process is responsible for the high 6 Li abundance 
in MPHSs, this is the limit to be adopted. Although pre- 
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vious studies 0,|2(| gi ye the 6 Li abundance of MPHSs as 
the upper limit, we adopt this as the lower limit consid- 
ering the possibility of stellar destruction. 

In Fig. |3| the contour of the lower limit l|29ll is shown 
by a solid line below the CMB constraint. Hence, a 6 Li- 
producing allowed parameter region certainly exists for 
t x = 10 8 - 10 12 s and Cx ~ 10~ 13 - 1CT 12 GeV. The 
parameter region allowed by the above constraints which 
also produces abundant 6 Li is marked as " 6 Li" . 

This 6 Li-producing zone is somewhat close to the limit 
deduced by fits to the CMB distortion. Thus, future 
measurements of the CMB spectrum may reach this pa- 
rameter region and provide better limits to this radia- 
tive decaying particle scenario. The Absolute Radiome- 
ter for Cosmology, Astrophysics, and Diffuse Emission 
(ARCADE) will observe the CMB at centimeter 

wavelengths with better sensitivity than current data. 
ARCADE may make it possible to impose a limit on the 
chemical potential up to \/j,\ < 2 x 1CP 5 and Compton y 
parameter up to y < 10 -6 , respectively. This sensitivity 
corresponds to the solid line marked as "ARCADE" in 
Fig- El It spans the parameter region for producing three 
times as much 6 Li as MPHSs (a dashed line, see dis- 
cussion below). Thus, if ARCADE does not detect the 
spectral distortion, 6 Li production from radiative decay 
can be ruled out. 

The allowed parameter region has a very small entropy 
increase of less than the order of 10 -5 , so that this has 
a negligible effect on the deduced rj. We have analyzed 
this parameter region to see the possibility of realization. 
It is important to check the results of the non-thermal 
nucleosynthesis for decay in the 6 Li-producing region. 
We made a model calculation with input parameters of 
t x = 1 x 10 10 s, Cx = 3 x 1CT 13 GeV and r) = 6.1 x 10" 10 . 
The final abundances obtained in this model are 



D/H 
3 He/H 
Y 

6 Li/H 
7 Li/H 



2.63 x 10"° 
2.48 x 1CT 5 
0.247 



= 4.69 x 10" 
= 4.36 x 10" 



-10 



(30) 
(31) 
(32) 
(33) 
(34) 



These are certainly consistent with the constraints we 
adopted in Sec. IIII Al The abundances of 3 He and 6 Li 
with respect to the SBBN abundances increase. In fact, 
the non-thermal production of 6 Li is unavoidably accom- 
panied by 3 He production. We confirmed that the pro- 
duced amounts of 3 He and 6 Li are proportional to (x 
when tx is fixed. However, for a longer decay lifetime, 
the 6 Li production is relatively more effective than the 
3 He production because of the different dependences of 
the production rates on the CBR temperature. 

If the inconsistency between the 7 Li abundance pre- 
dicted by SBBN and that measured from MPHSs is 
caused by depletion, 6 Li would have existed in the pri- 
mordial gas at a level larger than the abundance observed 
in MPHSs by at least the ratio of the SBBN 7 Li/H predic- 
tion to the mean value observed in MPHSs. The observed 




FIG. 4: Ratio of calculated 6 Li/H abundances (after the 
non-thermal nucleosynthesis) to the observed abundance in 
MPHSs as a function of tx- Results in the allowed parame- 
ter region of (tx, Cx) producing 6 Li/H larger than the value 
found in MPHSs, or the marked region " 6 Li" in Fig. [3] are 
plotted. The horizontal line indicates a factor of three over- 
production of 6 Li relative to the observed MPHS value of 
6 Li/H= 3 x 6.6 x 1(F 12 . The large circles denote values in the 
allowed region with abundances of 3 He/H=1.3 — 2.5 x 10 -5 
and 6 Li/H> 3 x 6.6 x 10~ 12 . The other parameters sets are 
indicated by small squares. 



7 Li/H abundance H| is 7 Li/H~ 1.62 x 10~ 10 . Hence, 
this factor is about (4.36 x 10- 10 )/(1.62 x lO" 10 ) - 3. 
So 6 Li should have been originally produced at an abun- 
dance more than about 3 times the presently observed 
value. We stress that the non-thermal 6 Li production 
inevitably brings about the production of 3 He, and this 
gives a strong constraint on the possible parameter space 
of unstable particles [l2L |20| . 

We have analyzed the upper limit to the 6 Li abundance 
resulting from the radiative decay process under the re- 
quirement of consistency with the other light-element 
abundances. In Fig. 0] the 6 Li abundances are plotted 
as a function of Tx- Points on this figure are allowed 
by the constraints imposed above and lead to 6 Li abun- 
dances above the level observed in MPHSs. The vertical 
scale is 6 Li/H normalized to the mean 6 Li/H abundance 
in MPHSs ( 6 Li/H)MPHS- The horizontal line indicates 
a factor of three enhancement in 6 Li. The large circles 
are for cases with more than three times as abundant 6 Li 
as the level found in MPHSs. Here, we adopt the one 
sigma 3 He/H=(1.9 ± 0.6) x 10" 5 [3l| as an extra con- 
straint. We note that, in cases adopting a tighter con- 
straint 3 He/H< (1.6 ± 0.3) x 10" 5 [Ij], one can still find 
an allowed region of tx = 3 x 10 10 — 3 x 10 11 s which sat- 
isfies the same constraint imposed on the 6 Li abundance. 
The small squares are for other cases of Eq. (I13|) . 

This figure confirms that 6 Li/H abundances as 
large as those in MPHSs multiplied by the ratio 
( 7 Li/H)/( 7 Li/H)MPHS can be produced by non-thermal 
nucleosynthesis without significantly impacting the other 
nuclide abundances. Although this explanation could re- 
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solve the discrepancy between the SBBN predicted 6 Li 
abundances and those derived from observations, it can- 
not resolve the lithium problem. This scenario neces- 
sarily requires some model for the stellar depletion of 
6 Li and 7 Li. Indeed, as discussed in [28j and references 
therein, models exist which suggest a very large deple- 
tion factor of 6 Li along with some 7 Li depletion. The 
production of 6 Li by radiative decay cannot explain the 
observed abundances of both 6 Li and 7 Li, if the stellar 
depletion proceeds as described by that model. How- 
ever, approximately equal amounts of depletion for both 
lithium isotopes can explain the measured abundances 
when combined with the non-thermal production of 6 Li. 
As for the case including the hadronic decay process 0] , 
it has been found that such particle decay could simul- 
taneously solve both the 6 Li and Li problem, even if a 
possible degree of depletion is included. Clearly, further 
research on such depletion processes, along with analysis 
of more observational data, is desired. 



VI. CONCLUSION 

We have investigated the possibility that non-thermal 
nucleosynthesis induced by the radiative decay of long- 
lived particles contributed to the light element abun- 
dances in the early universe. We have constrained this 
process by observed primordial light element abundances. 
The primordial nucleosynthesis induced by high energy 
non-thermal photons from the radiative decay of long- 
lived relic particles was calculated taking into account 
both the primary nuclear production reactions and the ef- 
fects of secondary production and destruction processes. 
We find, however, that the secondary destruction pro- 
cesses of primary D, T, 3 He and 6 Li have little influence 
on the final light clement abundances and in particular, 
it is confirmed that the effect of the destruction of the 
produced 6 Li is very small in the cosmic temperature 
regime where the non-thermal secondary 6 Li production 
is operative. 

We have utilized the observed light element abun- 
dances from various sources to explore the constraints 
on the parameter space of the lifetime (tx) and initial 
abundance (Cx) of unstable particles in a model in which 
the baryon-to-photon ratio r\ is fixed to the value inferred 
from the WMAP CMB power spectrum. We require that 
the non-thermal processes do not cause significant devia- 
tions from the observationally inferred primordial abun- 
dances. For each given lifetime Tx we deduce an upper 
limit to the product Cx of the photon energy of the radia- 
tive decay and the fraction of the number density of the 
decaying particles to that of the CBR photons. For short 
lifetimes, tx < 10 6 s, the lower limit to Cx is fixed by 
a D under-abundance, whereas for longer lifetimes, 3 He 
overproduction gives the strongest upper limit to Cx ■ 

The parameter values are also constrained by the in- 
duced distortion of the CMB black-body spectrum from 
energetic photons emitted in the decay. The CMB con- 



straint is more important than the light element abun- 
dance constraints for the case of a long lifetime tx > 
10 10 s. Future missions such as ARCADE would, per- 
haps, detect the signal of the radiative decay. Otherwise, 
the production of 6 Li in this process is constrained. 

Next we considered the free-streaming scale of the de- 
cay products. We imposed a constraint on the abundance 
parameter Cx under the assumption that the decay prod- 
ucts comprise a dominant constituent of the cosmological 
dark matter. For the specific model of Ref. [5^ a rela- 
tively late decay at tx > 10 5 s is forbidden. 

We studied the possibility of 6 Li production in non- 
thermal nucleosynthesis at a level which is compara- 
ble with the observed abundance in MPHSs. We find 
that there exists a parameter region leading to final 
abundances which are in reasonable agreement with the 
MPHS observations of 6 Li and is also consistent with 
the observational constraints on the other light nuclides. 
We show that if t x ~ 10 s - 10 12 s and Cx ~ 10" 13 - 
10 -12 GeV, it is possible that the radiative decay of long- 
lived relic particles causes the observed enhanced abun- 
dance of 6 Li in MPHSs. 

Analyzing the results of the light element abun- 
dances for the interesting particle parameter region in 
more detail, two important characteristics were found: 
One is that the excess 6 Li abundances from the non- 
thermal processes is regulated by the amount of 3 He co- 
production because 3 He is the seed for 6 Li in the pro- 
cess 4 He( 3 He,p) 6 Li. Hence, the radiative decay model 
which results in 6 Li-production above the MPHS abun- 
dance level is also reflected by an enhancement of the 
3 He abundance with respect to the SBBN value. There- 
fore, tighter constraints on the primordial 3 He abundance 
might exclude some parameter regions. However, a cer- 
tain range of the possible parameter region still remains. 

Another feature regarding this scenario of non-thermal 
6 Li production triggered by the radiative decay is that it 
does not resolve the lithium problem. Other mechanisms 
such as the stellar depletion of the lithium isotopes must 
operate to lower the 7 Li abundance in the atmosphere 
of MPHSs. In such mechanisms, however, 6 Li should 
be simultaneously depleted. In the parameter region, 
T X ~ 10 8 - 10 12 s and Cx ~ 10~ 13 - 10~ 12 GeV, the 
7 Li abundance is almost unchanged from the SBBN value 
even should the radiative decay of unstable particles and 
subsequent non-thermal nucleosynthesis be taken into 
consideration. Only 6 Li is enhanced by about three or- 
ders of magnitude higher than the SBBN value. This en- 
hancement can be as much as several times the observed 
6 Li abundance in MPHSs. Hence, the present scenario re- 
mains a viable possibility for 6 Li enhancement if the true 
depletion mechanisms or other systematic effects lead to 
approximately the same degree of depletion for both 6 Li 
and 7 Li. 

In summary, we have found a parameter region of 
tx ~ 10 8 - 10 12 s and Cx ~ 10~ 13 - 10~ 12 GeV where 
the non-thermal nucleosynthesis of 6 Li can explain the 
observed abundance level in MPHSs. This parameter 
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region satisfies the two observational constraints on the 
CMB energy spectrum and the primordial light element 
abundances, although it cannot be a solution to the small 
scale structure problem. 

It is very fascinating to consider the possibility that 
such new physics as the decay of exotic particles plays 
a role in producing the observed abundances of lithium 
in the universe. Clearly, further research into lithium 
nucleosynthesis and the interpretation of the 6 Li plateau 
would be very valuable. Also, further observations of the 
6 Li abundances in the stellar atmosphere of MPHSs will 
help to constrain the various suggested candidates for 6 Li 
production and to answer the question as to whether it is 
synthesized cosmologically or by cosmic ray interactions 
in the Galaxy. 
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